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Abstract 

Analytic solutions of the Regge- Wheeler equation are presented in the form of series 
of hypergeometric functions and Coulomb wave functions which have different regions 
of convergence. Relations between these solutions are established. The series solutions 
are given as the Post-Minkowskian expansion with respect to a parameter e = 2Mw, 
M being the mass of black hole. This expansion corresponds to the post-Newtonian 
expansion when they are applied to the gravitational radiation from a particle in circular 
orbit around a black hole. These solutions can also be useful for numerical computations. 
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1 Introduction 



In our previous workjQ], we presented analytic solutions of Regge- Wheeler (RW) equation 
in the form of series of hypergeometric functions. We proved that recurrence relations 
among hypergeometric functions as given in Appendix A in this text and showed that 
coefficients of series are systematically determined in the power series of e = 2Mw, 
M being the mass of black hole. We also presented analytic solutions in the form of 
series of Coulomb wave functions which turn out to be the same as those given by 
Leaver 0. We found that the series of solutions are characterized by the renormalized 
angular momentum which turns out to be identical. Then, we obtained a good solution 
by matching these two types of solutions. 

This method can be extended for Teukolsky equation|| in Kerr geometry. In this case, 
the coefficients of series of hypergeometric functions and also those of series of Coulomb 
wave functions series satisfy the three term recurrence relations. About these recurrence 
relations, OtchikfiJ] made an important observation that the recurrence relation for both 
series are identical, which enabled to relate these two solutions[|. Following the discussion 
by OtchikQ, Mano, Suzuki and Takasugi[||] extended our analysis to Teukolsky equa- 
tion in Kerr geometry and reported analytic solutions. We discussed the convergence 
regions of these series and the relation between two solutions of difference regions of 
convergences. The series are expressed in the e expansion which corresponds to the Post- 
Minkowskian expansion and also to the post-Newtonian expansion when they are applied 
to the gravitational radiation from a particle in circular orbit around a black hole. 

In this paper, we present analytic solutions of RW equation and discuss the analytic 
properties of these solutions by reorganizing our previous work|[| following our work 
on Teukolsky equation ||. Solutions are given in the form of series of hypergeometric 
functions and Coulomb wave functions and the convergence of these series are studied. 

1 In his paper, the relation between the series of hypergeometric functions and the series of Coulomb 
wave functions are related in the intermediate region where both series converge, though the series which 
he treated are not the solutions of Teukolsky equation. 
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The relation of these solutions with different regions of convergence are discussed. The 
e expansions of solutions are given up to the second order. It is not difficult to obtain 
solutions in the expansion of e which are used to estimate physical quantities in the 
gravitational wave astrophysics as discussed by Poisson and Sasaki 0. It is interesting 
to compare solutions given in e expansion and numerical solutions to test the accuracy 
of the e expansion. Since the convergence of series is fast, solutions can also be used for 
numerical computations. 

The Regge- Wheeler equation is given hjM 
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(1.1) 



where e = 2Mw, z = cur, M being the mass of Schwarzschild black hole and e the angular 
frequency. In the following, we summarize our result. 

In Sec. II, we derive analytic solutions of RW equation in the form of series of hyperge- 
ometric functions. They are X" n and X v out which satisfy the incoming boundary condition 
and the outgoing boundary condition on the horizon. These solutions are labeled by the 
renormalized angular momentum v — I + 0(e) which is determined by requiring the con- 
vergence of the series. It can be shown that —u — 1 is also the renormalized angular 
momentum. Solutions X\ n and X v out are invariant under the exchange of v and —v — \. 
These solutions can be expressed as combinations of other type of independent solutions 
X" and Xq v ~ 1 . These series solutions are convergent except for infinity as shown in 
Appendix B. 

In Sec. Ill, solutions are given in the form of series of Coulomb wave functions. They 
are Xq and Xq v ~ 1 . They are shown to be convergent for z > e as shown in Appendix 
B. These solutions can be expressed by linear combination of solutions which satisfy the 
incoming and outgoing boundary conditions, X^ in and X^ out . Because the three term 
recurrence relation among coefficients is the same of the hypergeometric function series, 
the renormalized angular momentum for the Coulomb wave function series is the same 
as for the hypergeometric case. 



In Sec. IV, we show that X v c is proportional to Xq. This is due to the fact that both 
series are specified by the same renormalized angular momentum v. The same is true 
for and Xq v ~ x and Xq v ~ 1 . Now we have analytic solutions both of which is expressed in 
the form of series of hypergeometric functions which is convergent except infinity and is 
expressed in the form of Coulomb wave functions which is convergent for z > e. 

In Sec.V, we discuss the relation between solutions of RW equation and solutions of 
Teukolsky equation. For some solutions, the relations can be explicitly proved. 

In Sec. VI, the Post-Minkowski (e) expansion of solutions are given up to the second 
order of e. The renormalized angular momentum v is expressed as a form of e expansion. 
We present the second order result. 

Summary and discussions are given in Sec. VII and the derivation of the three term 
recurrence relations among hypergeometric functions and among Coulomb wave functions 
which are crucial to derive the three term recurrence relation among coefficients are given. 

2 Analytic solution in the form of series of hyperge- 
ometric functions 

The RW equation has two regular singularities at r = 0, 2M and an irregular singu- 
larity at r = oo. In order to obtain the solution in the form of series of hypergeometric 
functions, we have to deal with these regular singularities. In particular, we eliminate 
the terms of 1/z 2 and l/(z — e) 2 in Eq.(l.l) by taking the parameterization 

X\ n = e^i-xY^l - x)-'pUx\ (2-1) 
where x — 1 — r/2M = 1 — z/e and e = 2Mu. Then, the RW equation becomes 

x(l - x)pl" + [1 - 2te + 2(1 + ie)x]pl n ' + (I - 1 - ie)(l + 2 + ie) P » n 
= _ ie[2x {l- x ) p ^ + 2x P i ; n + (l + ie)f m l (2.2) 
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We obtain the solution of this equation in the form of series of hypergeometric functions. 
For this, we introduce a parameter v and seek the solution in the form 

oo 

VU*) = E <Pn + ,(*), (2.3) 
n=— oo 

where 

T(n + v - 1 - ie)V{-n -v-2-ie) 
r(l - 2ie) 

x F(n + v - 1 - ie, -n - v - 2 - ie; 1 - 2ie;x). (2.4) 



Pn+u{X) 



As an initial condition, we set that a v n = 1 for n = and for others so that v = I + 0(e). 

In order for the coefficients of series (2.3) to be a solved, it is essential that the 
coefficients a^'s satisfy the three term recurrence relation. For this, it is essential that 
terms such as x{l — x)p' n+u and xp n+u are expressed as linear combinations of p n +u+i, 
p n+u and p n+u -i- This turns out to be true and the recurrence relations and the derivation 
of these relations are given in Appendix A. 

By substituting the form in Eq.(2.3) into the radial RW equation (2.2) and using the 
recurrence relations in Eqs.(A.7) and (A. 9) in Appendix A, we find that p\ n becomes a 
solution if the following recurrence relation among the coefficients is satisfied: 



« + i + « + TX-i = 0, (2.5) 



where 



v ie(n + v — 1 +ie)(n + v — 1 - ie)(n + v + 1 - ie) 

an (n + z/+l)(2n + 2z/ + 3) ' ( ' 

K = (n + ,)(n + , + 1) - 1(1 + 1) + 2* + (w + g < + *l + 1} , (2.7) 

v ie(n + v + 2 + ie) (n + v + 2 - ie) (n + z/ + ie) 
ln " (n + z/)(2n + 2z/-l) ' { ' j 

By introducing the continued fractions 

Rn(") = L n (v) = (2.9) 

a n-l a n+l 



we find 

*M = -ftT^bM- M " ) = -«T^M- (2 ' 10) 

From these equations, we can evaluate the coefficients by taking the initial condition 
Oq = 1. The parameter i/ is determined by requiring that the coefficients obtained by 
using R n {v) agree with those by using L n {y). This is the condition of the convergence of 
the series. This condition gives the transcendental equation for v 

R n {v)Ln-i{v) = 1- (2-11) 

The solution v of this equation is called the renormalized angular momentum. By using 
this v, the series is shown to be convergent in all the complex plane of x except for x = oo 
as discussed in Appendix B. 

The normalization of X" n is given on the horizon (x = e) as 

Y u f-ie f nV T(n + v-l- ie)T(-n - u - 2 - ie) \ ieH -x) f0 19 n 

Now we can show that —v—\ satisfies Eq.(2.13) so that —v—\ is also the renormalized 
angular momentum|5|]. As for the recurrence relation (2.8), we find aZn 1 = In an< ^ 
7Z^ _1 = a v n so that aZ^T 1 satisfies the same recurrence relation as a v n does. Thus if we 
choose <2q = aZ u ~ x = 1, we have 

< = aZ n - 1 . (2.13) 



By using the formula 

" r(n + z/ + 3-*e) 

x F(—n — v — 2 — ie, —n — v + 2 — ze; — 2n — 2v 



\— x' 



T(n + u-l- ie)T(-2n - 2v - 1) _ . n _ i/+1+ ^ 
+ r(-n-!/ + 2-ie) 1 J 

x Fin + z/ - 1 - ie,n + v + 3 - ie; 2n + 2v + 2; — 3—), (2.14) 

1 — x 
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we can show 

Xl = XZ + Xt v -\ (2.15) 

where 

x (l-x) n F(-n-u-2-ie, -n - u + 2 - ie;-2n -2u; ——). (2.16) 

-L Ob 

We consider the solution which satisfies the outgoing boundary condition on the 
horizon. The outgoing solution is parametralizd by 

Km = e^i-xTil - xYYouf (2-17) 
Then p v out should satisfy the following equation, 

*(1 - *)V V J' + [1 + 2^e + 2(1 - ie)x\tfj + (I - 1 + ie)(l + 2 - ie)f out 
= - ie [ 2x {l - x yj + 2(1 - 2ie)xp^ + (1 + 5ie)p v out ], (2.18) 

Now we expand p^ ut as 

oo 

= E <Pn + ,(^), (2.19) 
n=— oo 

where 

r(n + u-l + ie)T(-n-u-2 + ie) 
Pn+v[x) ~ r(l + 2^e) 

x F(n + u- 1 + ie, -n - u - 2 + ie;l + 2ie; x). (2.20) 

Similarly to the solution satisfying the incoming boundary condition, we find that the 
above series becomes a solution if a v n satisfies the same three term recurrence relation as 
a v n does. Then, by choosing <5q = a^^ 1 = 1, we find 

< = <■ (2-21) 
By using Eq.(2.21), the outgoing solution is given by 

oo 

Kut = e«*- X \-xf{\-x)- X £ <Pn + u(x) 

n=—oo 

= A^ + A^X^- 1 , (2.22) 



where 

A v = $^±M. (2.23) 
sin ir(v — it) 

This relation also shows that Xq and X ~ u ~ 1 are independent solutions of Eq.(2.1). 

3 Analytic solutions in the form of series of Coulomb 
wave functions 

Analytic solution in the form of series of Coulomb wave functions are given by 
Leaver P). Here, we parameterization to remove the singularity at r = 2M. By using a 
variable z = uor = e(l — x), we take the following form 

*c = (l--y e f v {z). (3.1) 

Then, we find 

2 ?f'J+[ x ?+2ez-l(! + l)]f v 

= e[z{fl + /„) - (1 - 2ie)f u -l(3 + 2te + e 2 )f u - ef v ]. (3.2) 

z 

If we consider v to be v = I + 0(e), the right-hand side of Eq.(3.2) is the quantity of 
order e so that this equation is a suitable one to obtain the solution in the expansion of 
e. 

Here we shall obtain the exact solution by expanding f v {z) in terms of Coulomb 
functions with the renormalized angular momentum u, 

U = t ^ r ^ + ^ 1 '^ + Z/+1 "" ) ^ +TO (^ (3.3) 
where F n+U is the unnormalized Coulomb wave function, 

F n+U = e-* z (2z) n+ »z T( r n + - ± 1 ± \ e) $(n + u + 1 + ie, 2n + 2v + 2; 2iz). (3.4) 

1 [2n + 2v + 2) 

By substituting Eq.(3.3) into Eq.(3.2) and using the recurrence relations in Eqs.(A.18) 
and (A. 19) in Appendix A, we find that the coefficients b v n satisfy the same three term 



recurrence relation as a v n does. Now we choose b v = b n 1 = 1 and then 

K = <■ (3-5) 

The fact that the recurrence relation for the Coulomb expansion case is the same to the 
one for the hypergeometric case means that the renormalized angular momenta v are 
the same for both solutions. This is important for these two solutions to be related each 
other. Therefore, —v—\ is also the renormalized angular momentum. As a result, another 
independent solution is obtained by replacing v with —v — 1. Thus two independent 
solutions are Xq and X^ v ~ x and Xq is given by 

xi = (i-ir i ■" ^riiT'i^r^-^ M. (3.6) 

V Zj n t^oo r ( n + V + 1 + *e)r(n + V + 3 + l€) 

Another forms of these solutions suitable for examining the asymptotic behavior are 

X C = X Cin + X Couti (3-7) 



y— v— 1 y—v—1 I y— v— 1 

= ie^ +1) ^ + X» Cm - ie^X^, (3.? 
sm 7r(f — «ej 



where 



x 



V n r(n + y - 1 - ie)T(n + u + 1 - ie) 



n=—co 



x m(n + v+ 1 -ie,2n + 2v + 2;-2iz), (3.9) 



X£ in = e~ lz z u+1 ( 1 - - ) ^ 2 u e-* e e m{u+1) 



oo 

X 

n=— oo 



T( 77 -1-7/ — 1 — If] 

E ' V(n + i/ + 3 + ie fo ( ~ 2 * )W * (w + ^ + 1 + ^, 2n + 2z/ + 2; 2iz). 

(3.10) 



In the end of this section, we present another set of solutions as 

X» c = (l-'-T t ^ ^ + ^-l + ^r(n + , + l-ie) 
C V ^/ n fr' 00 r(n + z/ + l+ie)r(n + z/ + 3-ie) 71 v+ny h V ; 



8 



where F u+n is defined in Eq.(3.4). In this parameterization, X c and X c u 1 are in- 
dependent solutions. There are two sets of independent solutions, {X c , X^^ 1 } and 
{X c , Xc"- 1 }. In the following, we consider a set {Xq, X c u 1 } which are related directly 
to a set of hypergeometric type solutions {Xq, Xq"^ 1 } given in the proceeding section. 



4 The relation between two solutions 

First we notice that Rq and R c are solutions of Teukolsky equation. Second we see that 
if we expand these solutions in Laurent series of 1 — x — z/en, both solutions give the 
series with the same characteristic exponent at x — > oo. Thus, Rq must be proportional 
to R c , 



Xq = K V X C . 



(4.1) 



The constant factor K v is determined by comparing like terms of these series. We find 



T(l + r + v + ie)T(— 1 + r + v + ie)T(3 + r + v + ie) sin7r(z/ + ie) 

~ T(n + u-l + ie)T(n + r + 2u + l) v 
x 2^ — ; — -a 



n=r 



X 



E 



(n-r)\T(n + u + 3-ie) 

T(n + v — 1 — ie)r(n + v + 1 — ie)a^ 
r(n + i/ + 1 + ie)r(n + i/ + 3 + ie)(r - n)!r(n + r + 2z/ + 2) 



(4.2) 



where r is an arbitrary integer. 

By using these relations, R v in can be written by using the Coulomb expansion solutions 

as 



(K v X Cin + K_ v _ x X c v in ) + (K u X Cout + K_ u _iX c u out ) 



K v + ie 



sin7r(z/ — ie) 



v-l 



X Cm ^{K v -ie^^K^ x )X c 



Cout' 



(4.3) 



The asymptotic behavior at z — > oo is 



(4.4) 



9 



where A u out and A v in are amplitudes of the outgoing and incoming waves at infinity of the 
solution which satisfy the incoming boundary condition at the outer horizon. They are 
given by 



A v 



(4.5) 



and 



AL 



g _ fe2 _i_ ie 



K v i v+1 + K-v-ii-i) 



, sin7r(i/ + ie) 



simr \v — ie) 



x 2^ K % 



(4.6) 



n=— oo 



r(n + i/ + 1 — ie) 

One application of these amplitudes is to derive the absorption coefficients. The 
absorption coefficient T can be expressed in terms of A m and A out as follows; 

2 



r" = i- 



A v 

^out 



A v 



(4.7) 



The upgoing solution which satisfy the outgoing boundary condition at infinity is 
given by X v Cout . This solution is expressed in terms of Xq and X U ^ 1 as follows 



Y{y - 1 - ie)T(v + 1 - ie) 
T(u + l+ie)T{u + 3 + ie) 



A Cout 



smniu + ie) 



sm7r \v — ie) 



-1 v-v-l 



K v )- x X»-ie™{K_ v _ x )- x Xz 



x 



e 2i™ + sinvr(z/ + ?e) 



sin7r(z/ — ie) 



(4.8) 



5 Relation between solutions of Teukolsky and Regge- 
Wheeler equations 

Let us present some solutions of Teukolsky equation in Ref.l for spin s = — 2 in the 
Schwarzschild limit (a — > 0) : 
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The solution which satisfies the incoming boundary condition on the horizon is ex- 
pressed as a series of hypergeometric functions as 



DC 



iC(- 2 ) = e i€X (-x) 2 - i£ a" n T F(n + v + l-iT,-n-v-ir;3-2ie;x) (5.1) 

n=— oo 

and coefficients a% T are determined by the three term recurrence relation, 

a n ' a n ' +1 + /5 n a n ' + 7n > o ft Li = 0, (5.2) 

„ )T ze(n + i/ - 1 + z'e)(n + u - 1 - ie)(n + z/ + 1 + ie) . . 

" n (n + z/ + l)(2n + 2z/ + 3) ' ^ ' ^ 

„ jT ze(n + z/ + 2 + ie)(n + f + 2 - z'e)(n + v - ie) 

7n ' (n + z/)(2n + 2z/- 1) ' [ ' 

and most importantly /3% is the same as the one in Eq.(2.9). By comparing the recurrence 

relation in Eq.(2.5) and (5.4), we find with the condition with clq T = ciq U ~ 1,t = 1 that 

A solution valid at infinity is given by 

R v c { -2) = z \ l -z) E H) n 



n=— oo 



Z, 

Yin + i/ — 1 + ze)T(n + z/ — 1 — ie)T(n + i/ + 1 + ie) 



X r(n + i/ + 3 + ze)r(n + z/ + 3 - ze)Y(n + i/ + 1 - ie) ^^^^ 
where 



Fl +V = e-(2z)^ r ^ + ^+J+ 2 Z ) £) $(n + „ + 3 + ie, 2n + 2u + 2; Kz). 

(5.7) 



Another solution is -^(-2) • These solutions are given in the expansion around the origin 
and are different from the ones given in our previous paper [T[] which corresponds to 
solutions expanded around the horizon. 

Let us first consider the following relation between a solution R of the Teukolsky 
equation for s = —2 and a solution X of the RW equation, 

n 



where r* = r + 2Mln(r/2M — 1) and A = r 2 — 2Mr. We substitute the incoming solution 
X? n in Eq.(2.1) for X in Eq.(5.1) and find 

R = ^ft;|-*4)V«"'(-*>-V fc 



= ± e »(*-V(- x y-xj!L p » 



dx 2tin 



n 



±ie( x -i) ( ) 2-ie T(n + v+l-ie)T(-n-v-ie) 
1 J „^oo r(3-2^e) 
xF(n + v + 1 — ie, — n — z/ — «e; 3 — 2ie; x), 



(5.9) 



where p" n is defined in Eq.(2.4). By substituting Eq.(5.6) into (5.9), we find that R in 
Eq.(5.8) is proportional to R u in ^_ 2 )- We find 

( d \r 2 ( d \ vu lV er(i/ + l-ie)r(-i/-ie)„„ 

(jF + «") a + -J rX - = e ^ r(3- 2fe) - R <«-v (5 ' 10) 

We next consider the inverse relation 

X = — - %oj —— %oj — — . (5.11) 

c A \dr* J A \dr* ) r 2 y ' 

For this relation, we consider a solution of the Teukolsky equation given in Eq.(5.4). 
Now, we substitute Rc(-2) into Eq.(5.7) and find 



2 



c dz z — e \ z 

= - * - —) 2 e"(z - er[e- iz (z - e)"* z~ 2 i$ ( _ 2) ] 

^r^+l+je) ~ r(n + 1/- l + »e)r(n + »/-l-»£) 

c i>+l-*e) 1 J r(n + z/ + 3-^e)r(2n + 2z/ + 2) n 

d 2 

x ^e- 2j2 2- 1 - te (2z) n+ ^<l)(n + v + 3 + ie, 2n + 2z/ + 2; 2iz) 
dx z 

u T(v+l + ie) e ie ~ I> + i/ - 1 + ?e)I> + // + 1 - it) 
c T(v+l-ie) { z ] r(ra + i/ + l + ie)r(n + i/ + 3-ie) n n+ "' 

(5.12) 

where F n+U is the Coulomb wave function defined in Eq.(3.4). The third equality is 
obtained by substituting a% T in Eq.(5.4) and the final equality is derived by using the 
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following relation 

d 2 

^[e- x x c - a+1 $(a, c; x)] = (c-a)(c-a + ^e^x^'^ia - 2, c; x). (5.13) 
ax 2 

By comparing Eq.(5.16) and Eq.(3.21), we find X is equal to X v c up to the numerical 
factor. Thus, we obtain the relation 



r 



qJA \dr* ~ lU ) A \dr* ~ / r 2 = « T(v + 1 - ie) Xc ' {5M) 

Equations (5.12) and (5.20) show the relations between a solution of Teukolsky equa- 
tion with s = — 2 and a solution of the RW equation. 



6 Post-Minkowskian expansions of solutions 

In this section, we discuss how to derive the solution in the expansion of the small 
parameter e = 2Mu. In order to find the solution in Eqs.(2.1) and (4.3) up to some 
power of e, we have to calculate v and a v n to that order by using Eq.(2.12) and (2.13) 
with the condition (2.14) and = a^ -1 = 1. Other coefficients b u n can be calculated 
from a v n by using the formula (3.10). 

For a v n with n > 1, the equation for R n (u) is useful. Since «^,7^ ~ 0(e) and and 
(3 u n ~ n(n + 21 + 1) ~ 0(1), we find 

R n (u) ~ 0(e) (6.1) 

for all positive integer n. As a result with = 1, we find 

< ~ 0(e") for n > 1. (6.2) 

Before discussing the coefficients for n < 0, we derive the renormalized angular mo- 
mentum v up to 0(e 2 ). For this, it is convenient to use the constraint for n — 1, 
i?i(i/)L (^) = 1. We notice that R\{v) ~ O(e) so that L (z/) must behave as 0(l/e). 
which requires that (3q + 7qL„i(z/) ~ 0(e 2 ) because «q ~ 0(e). In order to obtain v 
up to 0(e), we need to know the information of (3q up to 0(e 2 ) where the second order 
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term of v involves. Thus, we need the information about R\(v), L_i(z/), and 7q up 
to 0(e). Here we assume that L^iy) ~ 0(e) whose validity will be discussed later. In 
this situation, Ri(v), ckq and 7g can be calculated immediately. By substituting 

these to the constraint equation Ri(v)L (v) = 1 to find 



v = 1 + ' 



21 + 1 

^3 



o , -4 , (/-l) 2 (/ + 3) 2 {l-2f{l + 2f 



1(1 + 1) (21 + l)(2/ + 2)(2/ + 3) (2/- l)(2/)(2/ + 1) 



e 2 



+0(e 3 ). (6.3) 

In general, we can show that v is a even function of e. We see that the continued 
fractions R n (v) and L n (u) are determined by and o£t£ +1 . Since these quantities are 
even functions, the parameters in the transcendental equation R n (v)L n _i(v) = 1 aside 
from z/ are even functions of e which concludes that v an even function of e. Thus, the 
next correction term in v enters in the 4-th order term of e. 

The fact that the correction term of v starts from the second order term of e simplifies 
the calculation of the coefficients up to 0(e 2 ). 

Now we discuss the coefficients for negative integer n which are derived by using 
the equation for L n (v). For large negative value of | n |, L n (v) ~ ie/2n. Most of the 
negative integer value of n, L n (v) ~ 0(e). There arise some exceptions for certain values 
of n because the denominator of cx v n vanishes at n — — I — 1 and also f3% vanishes at 
n = —21 — 1 in the zeroth order of e. There is a speciality of Regge- Wheeler solutions 
that a y _ l+1 ~ 0(e 3 ). Because of these, we find 

L-i+M ~ 0(e 3 ) , 

~ 0(1) , 

L-2i-i{v) ~ 0(l/e) , 

^n(^) ~ 0(e) for all others. (6.4) 

From these observations, we find 

< ~ 0(e |n| ), for - 1 > n > -1 + 2, 
a-* + i ~ 0(e' +1 ), 
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a% ~ ~ 0(e l+2 ) 

< ~ 0(e w+1 ), for -l-2>n>-2l, 

< ~ CKe 1 ™ 1 " 1 ), for -2/-l>n. (6.5) 

With the above order estimates, we see that how many terms should be needed to cal- 
culate the coefficients with the specified accuracy of e. 

Comming back to v, we assumed that L_ 2 (v) ~ 0(e) which is valid since I > 2. . 

The coefficients a v n and also h v n up to 0(e 2 ) are obtained explicitly by 

-^ + 3) 2 (/ + 3) 2 2 3A 

01 = 2(Z + l)(2Z + l) e+ 2(Z + l) 2 (2Z + l) € +Q(£) ' (6 - 6) 

a2 ~ 4(/+f)(2/ + f)(2/ + 3)2 e + ° (e) ' (6 ' 7) 

°- 1 " 2/(2/ +f) e ~ 2/^(2/ + f) £ + ° (e) ' (6 - 8) 

-(/-3) 2 (/-2) 2 2 ^ 
^ 4/(2/-lV(2/ + V +Q(£); (6 ' 9) 
By using these coefficients, we can evaluate the incoming and the outgoing amplitudes 
in infinity. From Eq.(4.2), we find by taking r = that K v ~ 0(e~' _1 ). On the other 
hand, the estimate of K_ v _\ needs some care. By taking into account of the singular 
behaviors of gamma functions and the fact that the deviation of v from / starts from the 
second order of e, we find that K- v -\ ~ 0(e 1 ^ 1 ). Thus we obtain 



-~0(e 2i ), (6.10) 



which concludes that K- v -\ term contributes at most the order e 4 ( / > 2). In the 
approximation up to 0(e 2 ), we can safely neglect K_ v _\ term. 

Thus, we get the simple expressions for the outgoing and the incoming amplitudes as 
follows; 
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and 



where 



(6.12) 



B" = (-i) u+1 2- 1+ie e- m/2 . 



(6.13) 



By substituting the coefficients, we can easily calculate the amplitudes up to the order 
e 2 . Since the explicit expressions are complicated, we present the amplitudes up to 0(e) 
explicitly. We find 



A v 



I -it 



— e 



T(v - 1 - ie)T(v + 1 - ie) 
T(u + 1 +ie)T{u + 3 + ie) 

(/-!)/(/ + 3) + (/-2)(/ + !)(/ + 2) 



— e 



i 2 (i + iy 



i 



4(2/ + 1) 



21(1 + 1)(2/ + 1) 
(I - 1)1(1 + 3)(Z + 4) (_Z - 3)(Z - 2)(Z + 1)(Z + 2) 



(Z + l)(2Z + 3) 2 



z(2Z - r 



(6.14) 



A „ = {B5)T(u+l + ie, A 



out' 



(6.15) 



B v T(v + l-ie) 

The above result shows that the absorption coefficient V in Eq.(4.7) is zero up to the 
order of e for Schwarzschild black hole. 



7 Summary and Remarks 

We presented analytical solutions of Regge- Wheeler equation in the form of series of hy- 
pergeometric functions and Coulomb wave functions. The series are characterized by the 
renormalized angular momentum which turns out to be identical for both series which 
enabled us to relate these two series in the intermediate region where both series con- 
verge. The relation between solutions of Regge- Wheeler equation and those of Teukolsky 
solutions. This shows the consistency between present solutions and solutions presented 
in our previous paper p. 
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Solutions which we found will be used for the gravitational wave astrophysics, the 
gravitational wave emitted from a Schwarzschild black hole and also emitted from the 
particle rotating along a circular orbit around a black hole. Since solutions are given 
in simpler forms than those of Teukolsky equation, these solutions can be obtained in 
the e expansion in much higher accuracy so that these solutions can be used to test 
the e expansion. Solutions will be useful for numerical computations because the series 
converges fast. We expect that our solutions will become a powerful weapon to construct 
the theoretical template towards LIGO and VIRGO projects. 
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Appendix A: Proof of the recurrence relations among 
hypergeometric functions and Coulomb wave func- 
tions 

(a) Proof of the recurrence relations among hypergeometric functions 
Let us define 



A L , n = 



r(n + L - 1 - ie)Y(n - L - 2 - ie) 
Y{n + 1 - 2ie)n\ : 



(A.l) 



then 



We parameterize Pl+i as 



Pl+i 



Pl = J2 A L, n x r 

n=0 



™T(n + L- ie)T{n - L - 3 - ie) n 
2_ — : — ; ■ x 



n=0 

oo 



= E 



n=0 



1 + 



T(n + 1 - 2ie)n\ 
2(L + 1) 



n — L — 3 — ie 



so that we find 



E 



J n- L-3-ie 2(L + 1) 



(Pl+i-Pl)- 



Similarly, we find 



(A.2) 



(A.3) 



(A.4) 



E 

n=0 

Now we rewrite xpL as 



n + L — 2 — ie 2L 



^{Pl-Pl-i)- 



^ T{n + L-l-ie)T{n-L-2-ie) x 

XpL — / t ^ i - „ . s . X 



n=0 

oo 

= E 

n=0 

oo 

= E 

n=0 



T(n + 1 - 2ie)n\ 
n(n — 2ie) 



(n — L — 3 — ie) (n + L — 2 — ie) 



~rrAi, n x r ' 



1 + 



[L + 3 + ie){L + 3 - ie) (L — 2 + ie){L - 2 - ie) 



(2L + l)(n - L - 3 - ie) (2L + l)(n + L - 2 - ie) 



(A.5) 



(A.6) 
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By using Eqs.(A.4) and (A. 5), we obtain 



(n + v + 3 + ie) (n + v + 3 — ie) 



1 

+ 2 



2(n + z/ + l)(2n + 2z/+l) 
4 + e 2 



-Pn+u+l 



(n + v){n + v + 1) 



(n + i/ - 2 + ie)(n + i/ - 2 - ie) 
+ 2(n + z/)(2n + 2z/ + l) Pn+ ^ 1 ' 

Similarly, we rewrite x(l — as 



(A.7) 



x(l"lK = J2 nA L,n(x n -X n+1 ) 
n=0 

oo 

= J2 n 

n=0 

oo 

= 

n=0 

oo 

= -E 

n=0 



n — l)(n — 2ie) 



n — L — 3 — ie)(n + L — 2 — ie) 



(L + 3 -ie)(L + 2 + ie) (L - 2 + ie)(L - 1 - ie) 



4 + 



(2L + l)(n- L-3-ie) (2L + l)(n + L 
(L + 3 + ie)(L + 3 - ie)(L + 2 + ie) 



- 2 - ie) 



Al,nX n 



2L + 1 n- L 

(L-2 + ie)(L-2-ie)(L-l-ie) 1 



2L + 1 



By using Eqs.(A.4) and (A. 5), we find 



n + L — 2 — ie 



- 3 — ie 



(A.8) 



x(l -x)p' r 



(n + v + 3 + ie) (n + v + 3 — ir) (n + z/ + 2 + ie' 



n+u 



+ 



-2 + ie + 



2(n + z/+ l)(2n + 2z/+ 1) 

(4 + e 2 )(l + ?e) 
(n + z/)(n + i/ + 1) 



Pn+v+l 



Vn+v 



(n + v — 2 + ie) (n + v — 2 — ie) (n + v — 1 — ie) , . 

+ " J , w ^ 9 TTn -iW-i- A.9 

2(n + z/)(2n + 2v + 1) 

(b) Proof of the recurrence relations among Coulomb wave functions 



Let us define 



Br = 



"(2z) L z, 



Cl,u — 



T(n + L + 1 + ie) 
T(n + 2L + 2)n! ' 



(A.10) 
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then Fl is written by 



F L = B L "£C L , n (2iz) n . 

n=0 

Now let us consider the decomposition of Fl+i as 

Fh+l ~ 5i 75r(„ + 2L + 3)( n -i)! (2u) 



1 oo 

= *4e 



? " n + 2L + 2 



F L -2(L + 1)5 L ^ 



Ci, n (2iz) n 

C Ll n(2^) r 



n=0 



ra + 2L + 2 



Thus we have 



C L , n (2iz) 



^ Q n + 2L + 2 2(L + 1) 



(F L -iF L+1 ] 



Similarly, 



fi-i 



n=0 



r(n + 2L)(n)! 



= 5 L 2 



+ 1 ~ n + 2L + l 



n=0 



n + 1 



Then we find 



n=0 



T(L + ^e) 1 , C L>n (2iz) n 



+ 



_r(2L + l)2^ n + 1 
Now we evaluate Fl/z as 



— = %B L 

z 



2iBj 



^ r(n + 2L + 2)n\ W 
T(L + l + ie) 1 



= 2^ L 



T(2L + 2) 2iz ^ (n + 2L + 2)(n + l) 
f L + ie 



+ 



2L + 1 
L + l-ie 01 



-L „ — n 



2L + 1 ^ n + 2L + 2 



r(L + ^6) 1 ~ C L>n (2iz) n 
r(2L + l)2iz to 

' C L>n (2iz) 



n+1 
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By using Eqs.(A.13) and (A. 14), we find 



lp (n + u + l-ie) e 

z n+v (n + v + l)(2n + 2v + l) n+u+1 (n + v)(n + v + 1) " 



+ , + " + Wl , (A.17) 

(n + z/)(2n + 2z/ + l) n+I ^ 1 V ; 



Similarly, 



= -tfi + — *i + 2zS t £ + (A.18) 

Z n= Q Tl -+- ALi -+- z 

By using Eqs.(A.13) and (A. 14), we obtain 

P , _ (n + v)(n + v+l-ie) e 

(n + z/ + lj(2n + 2z/ + 1) [n + v)[n + v + 1) 

Appendix B: The convergence region of solutions 

From Eq.(2.13), we find 

li m n ^L = _ i im n ^L = _!£. ( B .l) 

By combining the large n behavior of hypergeometric functions^], we find 

lim w <P"+"( g ) = - l im n <Pn+Ax) ie _ ^ + _ 2 _ 1/2 
<_ 1 p n+I/ _i(a;) < +1 p n+I/+ i(:r) 2 

Thus the hypergeometric series converges in all the complex plane of x except for x = oo. 
As for the convergence of Coulomb series in Eq.(3.3), we findU 

lim = lim 1<*»M , = 2 -, (B.3) 



that 



l im < F "+»( Z ) = i im < F n+u{z) = _.e (B 4) 

<_!F n+l/ _i(z) n^-co a i +1 F n+u+1 {z) z 



Thus, the series converges for z > e or | x |> 1. 
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